Thermalization of coupled atom-light states in the presence of optical collisions 
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The interaction of a two-level atomic ensemble with a quantized single-mode electromagnetic field 
in the presence of optical collisions is investigated both theoretically and experimentally. The main 
focus is on achieving thermal equilibrium for coupled atom- light states (in particular dressed states) . 
We propose a model of atomic dressed-state thermalization that accounts for the evolution of the 
pseudo-spin Bloch vector components and characterize the essential role of the spontaneous emission 
rate in the thermalization process. Our model shows that the time of thermalization of the coupled 
atom-light states depends strictly on the ratio of the detuning to the resonant Rabi frequency. 
The predicted time of thermalization is in the nanosecond domain at full optical power and about 
10 times shorter than the natural lifetime in our experiment. Experimentally we investigate the 
interaction of the optical field with rubidium atoms in an ultrahigh-pressure buffer gas cell under 
the conditions of large atom-field detuning comparable to the thermal energy in frequency units. In 
particular, an observed asymmetry of the saturated lineshape is interpreted as evidence of thermal 
equilibrium of coupled atom-light states. 

PACS numbers: 42.50.Nn,05.30. Jp,32.70. Jz 



I. INTRODUCTION 

Present remarkable achievements with coherent ma- 
nipulation of coupled matter-field states evoke great in- 
terest in the investigation of phase transitions in such 
systems; see, for example, The key role in the 

behavior of coupled states under consideration is played 
by so-called dark and bright polaritons, that is, bosonic 
quasiparticles representing a linear superposition of pho- 
tons in an external (probe) field and the macroscopic 
(coherent) polarization of a two-level atomic system or 
excitons localized in quantum wells. 

The critical temperature of a phase transition for po- 
laritons can be high enough due to their small effective 
mass, which is many orders of magnitude smaller than 
the free mass of atoms (or electrons). Although evi- 
dence of Bose-Einstein condensation (BEC) of polaritons 
in semiconductor microstructures has recently been re- 
ported by several groups (see [l[ , , and Q ) observation 
of the high-temperature phase transition remains an un- 
solved problem. In this sense atomic systems seem to be 
more attractive and experimentally feasible for polariton 
BEC purposes [11]. 

The main difficulty with polariton condensate observa- 
tion is connected with the problem of achieving true ther- 
mal equilibrium for coupled matter-field states, which is a 
primary step in studying the phase transitions in the sys- 
tems under consideration. Roughly speaking, the polari- 
tonic system for the current experiments with semicon- 
ductor microcavities is in nonequilibrium (or quasiequi- 
librium) (see, e.g., @). One mechanism for achieving 
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thermal equilibrium of low-branch polaritons with the 
host lattice is to cool them with a phonon bath [Io| . 
The thermalization time must be shorter than the po- 
lariton lifetime in this case. For the current experiments 
the lifetime of polaritons is in the picosecond regime and 
comparable to the thermalization time. 

In the area of atomic physics, extremely long lifetimes 
of excitations are readily achieved at present. In this pa- 
per we show that the atomic polariton lifetime is limited 
by the lifetime of the two-level atomic transition T spont 
only, that is, by spontaneous emission. Since the value 
of T spont is in the nanosecond regime the thermalization 
time can be longer compared with that of polaritons 
formed in semiconductor devices. Therefore polaritons 
in an atomic physics system could be preferable for ob- 
servation of BEC, because of the long coherence times 
achieved. 

Progress toward the achievement of thermal equilib- 
rium in coupled atom-field states was made by some of 
us in In particular, the ability to thermalize coupled 
atom-light (dressed) states due to frequent collisions of 
rubidium atoms with buffer gas atoms in the presence 
of optical irradiation has been demonstrated experimen- 
tally. In the literature this process is called optical colli- 
sion (OC) (see, e.g., [HI)- We note that redistribution in 
OCs has also recently allowed laser cooling of ultradense 
atomic gases [l3j . 

In general atomic collisions in the presence of a laser 
field can be considered as a scattering (inelastic) process 
when both the internal (or kinetic) energy of the par- 
ticles and the energy of the scattered light are changed. 
Physically this leads to dephasing of the atomic polariza- 
tion that determines the broadening of the fluorescence 
spectrum and introduces an additional (collisional) phase 
shift. Notably, a closed representation to understand 
the effect of atomic collisions on the resulting width and 
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shape of spectral lines was performed by Weisskopf in his 
famous paper (l4j |. Thereafter collisional broadening has 
been a subject of intensive investigation both in theory 
and in experiment (see, e.g., |15l - [l8| ). Various aspects 
and approaches to line shape description, phase shift, 
and intensity due to nonresonant atomic collisions have 
been studied and summarized in [l9j . Notably, important 
properties of spectral line shift, width, and asymmetry 
for a wide range of foreign (buffer) atomic gas densities 
in connection with the so-called impact limit of atomic 
collisions have been discussed by A. Royer in [l8| (see 
also [l7|). 

A statistical approach to the OC problem was proposed 
in [20], based on the cross section of collisionally aided 
radiative excitation (or emission) for two-level atoms cou- 
pled to a thermostat of buffer gas. Subsequently, in [2lJ 
the authors established a simple theory of OC based 
on the dressed-state approach in the Schrodinger pic- 
ture. It has been shown that collisions with buffer gas 
particles reduce to a transfer between different dressed 
states of the atom. Spectral redistribution of quasireso- 
nant radiation occurring due to collisional relaxation has 
been a subject of experimental investigations as well (see, 
e.g., |22l - l24| ). The OC spectral line shape is essentially 
non-Lorentzian, that is, asymmetric. From a molecular 
point of view such an asymmetry can be understood by 
taking into account modifications of interaction potential 
curves [HHi]. 

Although the main features of OCs have been inves- 
tigated for a long time, the thermodynamic properties 
of coupled atom-light systems have not yet been studied 
fully. It it is pointed out in fl2j that, in the limit of a low 
Rabi splitting energy IiQr, that is, for hQji <C ksT (T is 
the temperature of the two- level atomic ensemble), the 
OCs reduce to equalizing dressed-state populations un- 
der the secular approximation. The deviations from the 
Einstein coefficients for absorption and stimulated emis- 
sion induced by OCs are discussed in (26|. However, the 
problem of thermalization of coupled atom-field states 
has not been studied. Standard theoretical approaches 
in this case are based on the rate equations for the pop- 
ulation of dressed states only, completely ignoring cou- 
pling between population and atomic coherences at the 
same time (see, e.g., [ll|, [Hj], and [11]). Such an ap- 
proach seems unsuitable if we keep in mind the problem 
of observation of polariton BEC. Actually spontaneous 
polarization buildup occurs in this case and polariton 
coherences become important (cf. @]). Thus, the role of 
the secular approximation in the thermalization process 
of coupled atom-light states should be clarified. 

The aim of this report is a fundamental theoretical and 
experimental investigation of the thermalization of cou- 
pled atom-field states, taking into account spontaneous 
emission processes beyond the approximations typically 
used. 

In Sec|TT]we establish the model of interaction of two- 
level atoms with the quantized optical field in the pres- 
ence of OCs. Realistic (experimentally accessible) condi- 



tions for OC of rubidium atoms with high-pressure buffer 
(argon or helium) gas atoms are discussed. In Sec lIIII thc 
Bloch-like equations for density matrix elements in the 
dressed-state basis are derived. Analysis of the steady- 
state solutions of these equations is performed in Sec lIVI 
We define the time of thermalization and specify neces- 
sary conditions to achieve thermal equilibrium for dressed 
atom-light states. In Sec|V] we derive the properties of 
intensities of spectral components for an atomic system 
in the presence of OC under the atom-field thermaliza- 
tion process. Thermalization of coupled atom-light states 
in an ultrahigh-pressure buffer gas environment is ex- 
perimentally investigated in Sec I VII The spectrum of 
fluorescence obtained from rubidium atoms is analyzed. 
We give an explanation for the observed resonance flu- 
orescence signal and necessary estimations of the time 
of thermalization according to the theoretical approach 
developed here. In Sec lVIII we summarize our results 
and discuss further prospects for the experiment and the 
corresponding theory. In particular, the possibility of 
polariton condensation is discussed. 

II. MASTER EQUATION APPROACH FOR 
ATOM-FIELD INTERACTION IN THE 
PRESENCE OF OC 

Atomic collisions in the presence of a nonresonant ra- 
diation field can be considered as the elementary process 
of a collision between an isolated two-level atom A and 
a foreign (buffer) gas atom of sort B; that is, 

A{a) + B + fwo L ♦=> A{b) + B, (1) 

which occurs with simultaneous emission (or absorption) 
of a photon with frequency uj l ■ The collision is called an 
OC when the frequency of the optical field ujl is near res- 
onant to an atomic transition in A. In FigfTJ schemes of 
nonresonant absorption of light with frequency u>l dur- 
ing collision with a buffer gas atom are shown. The 
case of collisionally aided excitation of level |b) of the 
atom for negative atom- light detuning, d — lol — luq, 
is illustrated in Figflja). For positive detuning 6 > 
(ujl > ojq), the atom has some energy excess HS, which is 
transferred to the kinetic energy of the atoms after the 
collision [FigOJb)]. Note that for both cases excitation 
of the upper level \b) is impossible due to atomic collision 
only. 

Our description of thermalization of atom-field 
(dressed) states developed here is based on the approach 
to OC presented in [12[ and generalizes their treatment 
to detuning of 5 values comparable to or higher than the 
thermal energy ksT . The master equation for the den- 
sity matrix cr in the presence of both OCs and radiative 
(spontaneous) relaxation processes can be written as fol- 
lows: 

da i . . fder] { do\ 
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resented as (cf. fl2j]): 
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FIG. 1: (Color online) Schematic representation of collision- 
ally aided absorption in a two- level atom with frequency u>l 
for (a) negative atom-field detuning, 5 < 0, and (b) positive 
atom-field detuning, that is, S > 0. The decay F is the spon- 
taneous emission rate. 



where the last two terms account for spontaneous emis- 
sion and collisions with buffer gas atoms. The Hamilto- 
nian H describes atom-field interaction under the rotat- 
ing wave approximation and has the form 



H = fkJ L pf + hu \b)(b\ + hg (S+f + SL/*) 



(3) 



where /(/') is the annihilation (creation) operator for 
the photons absorbed (or emitted) due to atomic colli- 



sions, g — \j \d a b\ uJL/2heoV is the atom-field interaction 
constant, which we take to be identical for all atoms, 
d a b is the atomic dipole matrix element, and V is the 
interaction volume. In expression Q S- = \ a )(b\ an d 
5+ = S_ = \b)(a\ represent atomic transition operators. 
The term 



dtr\ 

dt )rad 



-— (S+S-cr + aS+S-) 



(4) 



characterizes the contribution of spontaneous processes 
in Eq. @ ; T = 1/ T spont is the spontaneous emission rate 
(Tspont is the natural lifetime of the atomic transition). 
The last term in Eq.(|2]) describes the atomic collisions 
and can be established as 



^ Loll 
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cr + 2-fs z as z - ir\ [s z ,a] 



(5) 



where s z — | (\b)(b\ — |ffl)(ffl|) is the atomic population 
inversion operator obeying conditions s z \a) — — h\a) and 
s z \b) = ±\b). 

In ([5]) parameter 7 characterizes the collisional re- 
laxation rate (collisional broadening) in the presence of 
a monochromatic laser field; rj determines the average 
phase shift appearing due to collisions. In connection 
with the theory of OC parameters, 7 and 77 can be rep- 



7 = (1 — cos^ 
r\ = (sin < 



CO 



11 ' 

coll ' 



(6a) 
(6b) 



where <fi = J_°° [u>ba(t) — ojq] dt is the phase shift that 
accumulates during the collision. 

More rigorous expressions for collisional broadening 7 
can be found using a full quantum mechanical (micro- 
scopic) approach (see, e.g., |25|). In general, 7 implicitly 
depends on the atom-field detuning 5 and, thus, on the 
molecular potentials for a compound system. For 5 = 
the magnitude of 7 can be inferred from the expression 



7 ~ nplv T N B , 



(7) 



where po is the Weisskopf radius, depending on the 
molecular level variation (shifting) due to atomic colli- 
sions, vt = \f2ksT /m at is the thermal atomic velocity, 
and Nb is the number density of the buffer gas. 

In particular, in the experiment described in Sec lVIl 
we use rubidium atoms with mass m a t — 1-46 x I0~ 25 kg 
under the temperature T = 530K and a buffer gas of 
density Nb — 10 21 cto~ 3 . Taking into account typical 
values of the Weisskopf radius po ~ 10~ 3 pm for the col- 
lisional broadening rate ([7)1. one can obtain a value of 
few terahertz. The average collisional shift r\ can be ex- 
pressed via molecular level variation under the quantum 
mechanical description as well. Practically, r\ in (|6b[) has 
the same order of magnitude as a collisional broaden- 
ing 7. In the current experiment we used argon gas as 
the buffer gas (7/277 = 7^ r /27r ~ 7.2 GHz/bar, ri/2ir = 
riAr/2-rr ~ -6 GHz /bar) at 500 bar (50 MPa) pressure 
and helium as the buffer gas (7^ e /27r ~ 2.8 GHz/bar, 
ri He /2n ~ 1.6 GHz/bar) at 400 bar (40 MPa), respec- 
tively. 

The spontaneous emission rate for experimentally uti- 
lized rubidium atom D lines is T ~ 2ir ■ 6 MHz (T spont = 
27 ns) (see, e.g., [27j). Thus, condition 7 ~ 77 ^> T is held 
in the experiment under discussion. Nevertheless as we 
will see, spontaneous emission plays an essential role in 
the thermalization process of coupled atom-light states 
due to its ability to change the population of the excited 
atomic level. 

Note that the Doppler effect is not important for the 
problem under consideration. First, the motion of the 
atoms cannot change the population of atomic levels and, 
thus, cannot restrict the thermalization process. Second, 
Doppler broadening (of the order of a gigahertz) is sig- 
nificantly smaller than the collisional broadening in the 
experiment, thus neglecting the Doppler broadening in 
our calculations is justified. 

Let us now consider the situation when each collision 
happens in a short enough time span t co u. In this case 
two collisions are well separated in time; that is, we have 

Tcoll *C T co u, (8) 

where T co u is the time interval separating two collisions. 
Here we neglect temporal correlations between the spon- 
taneous emission process and the collisional one. Loosely 
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speaking spontaneous emission and OC optical collision 
processes are well separated in time as well. Further, 
equation ([2]) is valid under the so-called impact limit of 
OCs (see, e.g., 12]). The applicability of our model in 
this case is discussed later. 



III. BLOCH-LIKE EQUATIONS UNDER THE 
DRESSED-STATE REPRESENTATION 

In the absence of collisions with buffer gas atoms, 
Hamiltonian (J3J) for atom-light interaction explicitly has 
two eigenstates, called dressed states, defined as 

\l(N)) = sin9\a,N + 1) +cos9\b,N), (9a) 
|2(JV)> = cos0|a,iV + l) -sin.9\b,N), (9b) 

where N is the total photon number, and \a,N + 1) and 
1 6, N) are bare atom-light states. The mixing angle 9 = 
9(6) is defined by 



tan2<9 = < 29 < tt, 

o 



(10) 



where Oo ~ 2g\/N denotes the resonant Rabi frequency. 
In general, the parameters sin 9 and cos 9 can be repre- 
sented as [ll| 



sin a = 



cos 9 



± /i~X 

V2V Or' 



(11a) 
(lib) 



where flu — \J S 1 + f2p is the Rabi splitting frequency. 
It is important to emphasize that state \2(N)) is always 
located below state |1(JV)). 

In Fig[2]the dependence of dressed-atom-state energies 
as a function of light frequency ojl is presented schemat- 
ically. The frequency gap between dressed states | l(iV)) 
and 1 2 (TV)) is the Rabi splitting frequency The gap 
is minimal and equal to the Rabi frequency Qn :m i n = Q Q 
under the atom- light resonance condition (5 = 0). In 
this paper we mostly focus on the case for which the 
Rabi splitting energy is comparable to or higher than 
the thermal energy; that is, 



(12) 



For the current experiment (see Sec I VI I) we are limited 
to Rabi frequencies f2o/27r by the value of 0.1 THz, which 
corresponds to full optical power P ~ 300 mW. At the 
same time the thermal energy (ksT) for rubidium atoms 
at ambient temperatures (T = 530 K) corresponds to a 
frequency of 11 THz. In this case, we need to have large 
detuning 5 values (|<5|/27r > 11 THz) to fulfill relation 
(|12[) . In other words, in this paper we are practically 
interested in the perturbative limit when 




FIG. 2: (Color online) Dispersion relation for coupled atom- 
light states versus u)l', is the resonant Rabi splitting en- 
ergy. The upper branch corresponds to the |l(iV)) dressed 
state, the lower branch describes the \2(N)) state. The dashed 
horizontal line corresponds to the uncoupled level \b, N), that 
is, to the atomic transition energy !ko . The experimentally 
observed coupled states are states (2) and (3), for which con- 
dition (Tl3l) is satisfied. 



Physically condition (|13l) means that we deal with 
dressed states situated far from the region of resonant 
atom-field interaction. In Figf2] this situation is labeled 
as (2) and (3) respectively. 

We solve the master equation ([2} in the basis of dressed 
states |l(iV)) and \2(N)). The density matrix elements 
in the dressed-state representation © traced over the 
photon number N are defined as 



an = 53<l(JV)k|lW) 5 (14a) 

JV 

a 22 = J2(^NMHN)), (14b) 

JV 

«7 3 i=crI 2 =5^(2(JV)Hl(JV)). (14c) 



N 



The matrix elements an and 022 describe the popula- 
tions of the dressed states |1(AT)) and \2(N)), respec- 
tively. The nondiagonal elements o\2 (°2i) characterize 
dressed-state coherences and correspond to population 
transfer between dressed-state levels. Notably, the to- 
tal population of dressed states ([T4l is conserved for OC 
processes; that is, we have d (an + 022) / 'dt = 0. 

We now consider the properties of the real components 
S x ,y,z of the pseudospin (Bloch) vector S combined from 



matrix elements (1141) as follows 



S x = C12 + (721, 



S y = i (<7i2 - 021) , 

S z = <J\l — 022- 



(15a) 
(15b) 
(15c) 



Oo« \S\. 



(13) 



With the help of ©-© it is possible to get Bloch-like 
equations for pseudospin components S Xl y }Z , 
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dSx 
dt 



(eq) 



(Xcoh +$)S X - QrS v - 2 (a - Yyi) S z - 



dt 



dSy 

~dt 



= -2w 



where the following notations are introduced: 



2wS, 
(2w + r+) 



Tsin(20), 



(eg) 



1 + I S z - 



2(a-r 12 )S x -2USy + T-, 



(16a) 
(16b) 
(16c) 



7sin(46>) 



7 sin 2 (29) 



S = 7Cos 2 (26»), [/ = 



T co/l = - (1 + sin 2 (2#)) , r± = r (sin 4 (#) ± cos 4 (fl)) , r 12 = 



r/sin(20) 
2 ' 
T sin(40) 



(17a) 
(17b) 



r 



\ t 



i->i 



,r r 

1 A 2->2 x 2->ll 



l 22 



1— >2J 



I t 



|2(JV)> 



\1(N-1)} 
\2(N-l)} 



FIG. 3: (Color online) Energy (population) transfer between 
levels of the dressed atom induced by dephasing collisions. 7o, 
In and I22 indicate intensity weights of three-color photons 
emission. In addition, spontaneous emission described by the 
rates Fi_>j = 1,2) causes a transfer between levels of 

different manifolds. 



In Eqs. flTB")) . £Ir = fin — Srj/flu is a modified Rabi 
splitting that takes into account the average phase shift 
7] arising due to collisions with buffer gas atoms. The 
coefficients T co h, Y±, and Ti2 in (|17b[) characterize spon- 
taneous emission in the dressed-state representation. In 
general they are connected with spontaneous emission 
rates = T 2 ^ 2 = fr sin 2 (20), r 2 _>i = rsin 4 (6>) and 

Ti^2 = rcos 4 (#) from dressed-state levels linking neigh- 
bor manifolds (see Fig(3J). U describes the influence of 
the average phase shift on the evolution of the system 
under discussion. 



Since we are interested in the conditions for the tran- 
sition of coupled atom-light states to thermal equi- 
librium, we have introduced in Eqs. (|16l) the term 
2wS l f qS> / (2w + T + ) for the evolution of the Bloch vec- 
tor component S z . The presence of such a term in the 
set of Eqs. (fTrj)) is compatible with existing theories of 
spin-boson interaction with thermostat particles (due to 
collisions in our case) discussed by A. Leggett et al. in 
(28j . In particular, according to our approach in the ideal 
case, neglecting spontaneous emission terms in (|16|) . the 
dressed-state population imbalance S z should reach its 
thermodynamically equilibrium value S z eq ^ with the rate 
2w due to collisions with buffer gas atoms [see (l2"Tj) ]. 
which is in agreement with experimentally tested ap- 
proaches to OCs based on the solution of Boltzmann-like 
(rate) equations [Til. l26j . 

The dependence of Sz on temperature can be easily 
understood from the thermodynamic properties of the 
coupled atom-light system (cf. [ll|, [.29]). In particular, 
the population of the lower dressed state |2(7V)) should 
be much larger [exp (SQft/fegT) times] than that of the 
upper one in thermal equilibrium. We suppose 

= ~ tanh (hn R /2k B T) sa - tanh (ft \S\ /2k B T) , 

(18) 

where the latter expression is valid under the perturba- 
tive limit (|13[) . In particular, for a near-resonant atom- 
field interaction (when inequality h\S\ ~ hQn <C k B T 
is still true), the equilibrium value of S z according to 

Eq. (fl8|) is S z eq ^ — and we arrive at the well-known re- 
sult for which collisions with buffer gas atoms tend to 
equalize the population of the dressed states (cf.pjj]). 
In the opposite limit of large detuning one can omit 
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the terms containing collisional broadening 7 and spon- 
taneous emission T from Eqs. (ll6[) and obtain the same 
results as in [i^ 1 obtained from the Schrodinger represen- 
tation for OCs. In particular, without emission the colli- 
sions do not induce transitions between state \a, N + 1) 
and state \b,N). Collisionally aided excitation allows a 
transfer between dressed-state components also in the 
case of large detunings. The energy difference is balanced 
by the kinetic energy of the colliding particles. Since 
state \2(N)} is energetically lower than state |1(7V)), en- 
ergy is taken from the thermal reservoir of the buffer 
gas during the transition from |2(JV)) to \1(N)) (Fig©. 
About 10 3 -10 4 collisions happen during the natural life- 
time T S pont of rubidium atoms. 

The population transfer between dressed states evokes 
a thermalization process of coupled atom-light states 
which is characterized by the thermalization rate 2w = 
7 (Qg/Q^) fa 7 (^o/<5 2 ) in Eq.([T5c]l for the dressed-state 



where we made the denotation D = Q R + (T co h + 

0[(r/2) + 7 ]- 

First, we examine the role of atomic collisions in the 
thermalization process, completely neglecting sponta- 
neous emission within the so-called secular approxima- 
tion: 

r < 7, t) < n R . (20) 

In this case from expressions (|19[) one can assume that 

s x st) « 4 st) ~ °> ( 21a ) 

Si st) ^S ( z eq) . (21b) 

The inequalities (|20|) represent the necessary condition 
for achieving a thermodynamically true equilibrium (|21b ) 
for the dressed-states population. The coupling between 
dressed-state coherences <7i2 (°2i) and populations crn, 
CT22 can be completely neglected in this case; see (|21b ) 
(cf. 0). 

The thermodynamically full equilibrium behavior of 

the dressed-state population imbalance Sz is repre- 
sented by the dotted (red) curve in FigfJ] It is important 
to emphasize that the atomic ensemble under discussion 
indicates a two-level system without inversion under the 



population imbalance S z . Notably, 2w depends on the 
collisional rate 7 as well as on the ratio £l /6, which char- 
acterizes the atom-field interaction. The dependence of 
2w on the laser intensity is in agreement with present 
theories on OCs (see, e.g., [251]), whose validity in the 
yet unexplored high-pressure buffer gas regime with a 
high multiparticle collisional rate, however, remains to 
be tested. 



IV. THERMALIZATION OF COUPLED 
ATOM-LIGHT STATES 

Our goal in this section is to find stationary solutions of 
Eqs. pljf which are close to the thermodynamically true 
equilibrium state S z eq ^ of the coupled atom-field system. 

The full set of Eqs. ([T6|) yields steady-state solutions: 



(19a) 
(19b) 
(19c) 



perturbative limit (|13p for negative atom-field detuning 
6 only. Actually in this limit from dressed-state defini- 
tions @, we have sin 9 ~ 0, cos 6* ~ 1, and the lower 
dressed state \2(N)), which is macroscopically occupied, 
corresponds to the bare state \a,N+ 1), which describes 
atoms in the ground state. 

For large positive atom-light field detuning 6 > under 
condition (fT5|) , we can put sin 8 ~ 1 , cos 9 ~ , in Eqs . © , 
which implies that the lower dressed state \2{N)) corre- 
sponds to the excited atomic level \b, N), which is much 
more populated. In this case we achieve inversion in the 
two-level atomic system being under thermal equilibrium 

(cf. Hi). 

The decay of atomic upper level \b) with the rate T 
leads to nonequilibrium processes in the dressed-state 
population behavior. Taking into account (|16b ) spon- 
taneous emission under the assumption (f2Ta ). we get a 
steady-state solution for the dressed-state population im- 
balance S z : 

2 W (e)si eq) +r_(g) 
2w(0) + r + (6) ■ [ZZ) 

In FigH] the dependence of the dressed-state popula- 
tion imbalance S z in the steady state (|2"21 on the atom- 
field detuning 6 is plotted. The full thermalization of the 
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T S m2e(^+^ +2(f/Q R +(a-r 12 )(^ +7)) ^ 
2(U(T coh + <;)- Q R (a - T 12 )) I S< rt > 
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(eg) 



2wS 



(eg) 



s (st) _ 



2wS 



(eg) 



2w + r+l 

DT- - 2rsin20[(a - T u )(l + 7) + UQ a ] 



2w + r+ 1 

- n R T sin 29 
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atom-light state (dotted curve in Fig|4j is approached 
at infinite Rabi frequency (i.e. for infinite input laser 
power) only. The same result can be obtained from 
Boltzmann-like (rate) equations for dressed-state popu- 
lations if the contribution from spontaneous emission to 
the thermalization of the atom-field states is completely 
neglected (cf. (llj). The asymptotic behavior of popu- 
lation imbalance presented in FigU for large values of 
atom-field detuning S can be easily analyzed in the per- 
turbative limit (|13[) under condition (|12p . In particular, 
from (|2"2"j) one can obtain 

"n _ 2 e -W\/kBT 

S z w Si eq) + ^ — for S < 0, (23a) 



1 , -r n o 



S z « + 



1 + r<5 2 



for 5 > 0. (23b) 



The terms containing T in (|2"21 and (j2"3")l characterize 
the influence of spontaneous emission on the thermal- 
ization process. The decay of the dressed-state popula- 
tion imbalance in (|22[) due to spontaneous emission is 
described by T + . To minimize this effect we require the 
rate of thermalization of atom-field dressed state 2w to 
be much higher than the effective rate of spontaneous 
emission T + ; that is, the inequality 



2w > r. 



(24) 



should be satisfied. 

With the help of definitions (JTTJ) and (|2"4"1) . one can 
obtain 



r 



- < < 1 



7 



5- 



(25) 



Relation ([25|) is one of the main results of the paper 
and a characterizing condition for reaching true ther- 
mal equilibrium of coupled atom-field (dressed) states. 
Such an equilibration is achieved at the far detuned tails 
of dressed-states population imbalance S z in FigJH for 
which, with condition ([2"3")l . we have S z « S z eq \ 

The suppression of the thermalization process due to 
spontaneous emission leads to the formation of ther- 
modynamically quasiequilibrium of coupled atom-light 
states for which condition (l2"5l) is violated. Indeed, for 



nl r 1 

o z 7 



(26) 



from ([23|) we get 



Sz « + ^ - 2e- n l 5 l/ fcsT for S < 0, (27a) 

S z w S"^ 9 ) + 2 for 5 > 0. (27b) 

For highly nonequilibrium coupled atom-light states 
S z ~ 1 [see ([2Tb ) and FigH] and spontaneous emission 
decay drives the dressed-state system out of equilibrium 
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FIG. 4: (Color online) Dependence of dressed-state popula- 
tion imbalance S z as a function of detuning 5/2n for 500-bar 
argon gas collisional broadening y/2ir — 3.6 THz; collisional 
shift 7//27T = —3 THz. The resonant Rabi splitting frequency 
f2o/27r is infinite for the dotted (red) curve, 0.1 THz for the 
solid (black) curve, and 0.03 THz for the dashed (blue) curve. 
In all cases the gas temperature is T = 530 K and the spon- 
taneous emission rate is T ~ 37 MHz. 



for large positive detuning S. In conclusion, Eqs. pli 
yield a stationary solution only. 

Conditions (f2"4"|) and (f2"5)) can be represented in some 
other form introducing the time of thermalization T t herm 
of coupled atom-field states as follows: 



2tt 
2^ 



2tt 

7 



1 



O 2 

"0 



(28) 



According to this definition the minimal time of ther- 
malization y^ min ) 
atom-field interaction 



therm — ^ n / 1 occurs for a near-resonant 
-C fi ) representing state (1) 

in Fig. 2. The time r/^-m ^ s m tnc picosecond regime 
for experimentally accessible collisional broadening 7. 

For large atom-light detuning in the perturbative limit 
(|13p . the time of thermalization Ttherm increases and ap- 
proaches 



Tfhe 



2tt<5 2 2nk 2 B T 2 



(29) 



Notably, the quantity T t herm is inversely proportional 
to £Iq, that is, to the incident optical power Pq. Taking 



into account definition (|28[) . we can express condition 
for thermalization of atom-field coupled states as 



T, 



I lie 



spont • 



(30) 



Thereby, at appropriate time scales the thermalization 
time Ttherm of the atom-field dressed states must be 
faster than the atomic transition lifetime T spont . In par- 
ticular, for our experiment the minimal achievable time 
of thermalization at full optical power is about 3.37 ns 
for argon buffer gas and 10.8 ns for helium (see Sec lVIl) . 
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Now let us examine nonequilibrium properties of the 
atom-light states. The value of detuning 5 is of the order 
of the collisional broadening 7 in our experiment. In fact, 
we are working beyond the secular approximation limit 
established by (l20l) . In this case the Bloch components 
S x and S y can be adiabatically eliminated from Eqs. (|16p 
and we can obtain the following solution for the dressed- 
state population imbalance S z : 

S z (t) = Si st) + (5,(0) - e ^(2-+ r +) t , (31) 

where S z (0) is the initial value of the pseudospin compo- 
nent S g at t = 0, and Sg is determined by Eq.(flT)b). 

The numerical solution of the full set of Eqs.(fT!)|) re- 
vealing the quasiequilibrium dynamics of the coupled 
atom-light system for 500-bar argon buffer gas is pre- 
sented in FigJS] We suppose that all atoms initially oc- 
cupy the lower dressed-state level |2(7V)); that is, S z (t = 
0) = — 1. FigureEJa) demonstrates the decreasing atomic 
polarization described by dressed-state coherences, the 
Bloch vector components S x and S y , respectively. In par- 

(st) 

ticular, in the steady-state S x , y tend to the values S x , y 
according to Eqs. (fT§a ) and (|19b ). They clearly exhibit 
the presence of some small (residual) polarization of the 
atomic medium which depends on the vital parameter 
T + /w. The dressed-state coherences S x and S v are rela- 
tively large for nonequilibrium states when T + 3> w. On 
the contrary, we can use Eq. (|2Ta ) and ignore the dressed- 



state polarization at full thermal equilibrium under the 
secular approximation ([20)1 . 

The transient regime of thermalization of coupled 
atom-light states for the experiment described in Sec lVIl 
is presented in FigJSJb). The dependences based on a 
numerical solution of full Eqs. (|T6|) are represented by 
solid curves in FigJSJb). Dashed curves correspond to 
an expression (|3"Tj) . The steady-state levels for various 
resonant Rabi frequencies £1 /2tt are shown by horizon- 
tal (magenta) lines. These levels are determined by the 
magnitude of population imbalance S z for a detuning 
\5\ = ksT /H that corresponds to a frequency of about 11 
THz. From FigJSJb) it follows that for sufficiently long 
time scales such as t ^> l/(2w + r+), the dressed-state 
population imbalance S z (t) reaches its steady-state value 

S { z st) (for a given Rabi frequency Q /2ir) as determined 
by expression (|T5b ). 



V. INTENSITIES OF SPECTRUM 
COMPONENTS 

Here we discuss the modification of the intensities of 
the fluorescence spectrum components under the ther- 
malization of atomic dressed states. In general the flu- 
orescence triplet consists of a central line of frequency 
wl and two sidebands centered at wl ± flu (see FigJ3|). 
These lines have intensities with weights 



i n = i{uo L + n R ) = ffu r^ 2 = r cos 4 (0)- 



w(9) 



1 + S : 



(eq) 



r sin 4 e 



2w(0)+T+(0) 



I22 = I{^l - Or) = o-22r 2 ^i = rsin 4 (6 l ) 



w(0) 



(eq) 



T COS 4 9 



2w(0)+r + (6) 
I = I{lo l ) = T sin 2 {6) cos 2 (6). 



r 



(32a) 

(32b) 
(32c) 



Following (j32fc) the la intensity component weight is 
still unchanged during the thermalization process be- 
cause it does not depend on temperature. In other words 
thermalization of coupled atom-light states reduces to 
redistribution of intensities between the In and the J22 
components. 

To be more specific, we examine Eq. (|32l) in the per- 
turbative limit (fT3|) for negative detuning 5 < 0. In this 
case we can assume that sin 2 (#) w Hq/45 2 , cos 2 (0) ss 1 
[see (jlip]. At the same time in this limit we have 
|1(JV)) - \b,N) and \2(N - 1)) - \a,N) for dressed levels 
of one manifold. Physically this means that In approx- 
imately corresponds to the transition from the upper to 
the lower state (see Fig|3]). 

Let us now analyze Eqs. (|32[) in the limit where ther- 
malization does not occur, that is, when condition (|26[) 



is fulfilled for coupled atom-light states. From 
get 



" 16^4 ' 

Io - W 



we 



(33a) 
(33b) 
(33c) 



In expressions (I33[) we arrive at familiar results (see, e.g., 
[l2l |) in the limit of near-resonant atom- field interaction 
when h\S\ <C ksT and In sa 2/q7 /r Iq. However, in 
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FIG. 5: (Color online) Dependences of (a) dressed-state pseu- 
dospin components S x , y and (b) population imbalance S z as 
a function on reduced time yt. Parameters are: "//2tt — 3.6 
THz, r?/27r = -3 THz, T ~ 37 MHz, 8/2ir = -11 THz, and 
T = 530 K. The resonant Rabi splitting frequency flo/2-K is 
0.1 THz for curves (1) and 0.03 THz for curves (2). Initial 
conditions are: S x , y (0) = and S z (0) = — 1. The vertical line 
corresponds to the time Ttherm of thermalization of coupled 
atom-light states for Q.q/2-k = 0.1 THz. 



the general case, from (13"3")) we obtain 
hi Aye~ h ^ /kBT 

"io~~ f • 



(34) 



For large detuning 5, one can have In <C Iq under the 
condition K\6\ ^> /c^Tln [47/r]. However, such a con- 
dition is not realized for the experimentally accessible 
range of detuning 8 and temperatures T - cf. (fP2"j) . 

We can show that in the presence of thermal equilib- 
rium, under the conditions (|24|) and (|25|) the In intensity 
component changes essentially, and from (|32b ) we obtain 



.-(therm) 



Te -h\s\/k B T^ 



(35) 



At the same time the other intensity weights, I22 and Iq, 
still remain unchanged [see Eqs. (|3"3"b) and (|33b)]. From 



33k) and (|35| it is easy to see that 

ill < Jn ) 



(36) 



which implies a significant increase in In under the ther- 
malization process. 

Now let us turn our attention to positive-valued large 
atom-field detuning, that is, to 8 > . Without thermal- 
ization (2w < r), from ([32k) and (|32b) we get [cf. ([33k) 
and {33))]: 



/11 f= 
^22 



16(5 4 ' 
' <5 2 ' 



On the contrary, Eqs. ffiZh ) and (|32b) approach 



jitherrci) 
'll 



ro 4 

16£ 4 



22 



r, 



(37a) 
(37b) 

(38a) 
(38b) 



under the fulfillment of conditions (f2"4")l and (|25l) . 

Thus, thermalization of coupled atom-light states for 
positive detuning 8 is mostly due to dramatic increase 
in the I 2 2 intensity weight component and suppression of 
the In component simultaneously. 

It is also fruitful to introduce the total intensity weight 
of fluorescence (absorption) J as a sum of all intensity 
weights represented in (|32|) ; that is, I = In + I22 + h- 
With the help of definitions © and flU), from fl3JJ| we 
get 



I = a bb T = - 



r cos 2 (26») \ Twcos(26)S, 



2w 



2w + T 4 



where abb — (b, N\a\b, N) is the population of the upper 
(excited) atomic state. 

The thermalization of the atom-light field states is de- 
termined by last term in (l39l) (cf. [12j]). In Fig J6] the 
dependence of the reduced intensity weight J/F (excited 
atomic-state population a b b) is presented as a function of 
atom-light field detuning 8. The asymmetry of the be- 
havior of / in Fig|5]can be easily understood by introduc- 
ing the difference AI between positive (far blue wing)- 
and negative (far red wing)-valued detuning intensities 
for which, from (|39|). we obtain 

2Twta,nh(h\S\/2k B T) 



AI = I (8 > 0) - I (8 < 0) 



2w + r H 



(40) 

For complete thermalization of the atom-field states 
[dotted red curve in FigJB], under condition (|2"5j) the in- 
tensity weight of fluorescence / can be approximated by 
a Fermi-Dirac distribution function as 



(41) 
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FIG. 6: (Color online) The reduced intensity of fluorescence I 
versus atom- field detuning S/2it for 500-bar argon buffer gas 
at different values of the resonant Rabi frequency f2o/27r. 
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creases in the incident optical power, that is proportional 
to the square of resonant Rabi frequency SIq [cf.([3"9])]. 

To end the theoretical part of the paper, let us briefly 
discuss the boundaries of applicability of our theoretical 
model for OCs. It is necessary to accentuate the impact 
limit (QrT cM <C 1) of OCs when the processes of Rabi 
oscillations, spontaneous emission, and atomic collisions 
are mutually uncorrelated. Only in this case is the master 
equation ^ used in this paper valid. In our experiment 
the impact limit is not completely fulfilled (see SeclVip. 
In particular, only inequality 51 -C t~ o11 takes place in 
our case. At the same time the experimentally accessible 
atom-light detuning 6 is of the order of the inverse time 
of collisions t~1. 

Further, we do not take into account the influence 
of the Franck-Condon factor (overlapping integral) in 
electric dipole operators, which characterizes overlapping 
between the (quasi-) molecular wave functions for ini- 
tial and final states during the collision (see, e.g., j3(ij|). 
Obviously, this factor becomes more important when 
the impact limit breaks down (for detuning S such as 
|<5|t coJ; 3> 1) and a full quantum mechanical approach 
to OCs and line profile description become necessary (cf. 
[Hi, [3l|, [11|)- Nevertheless, as we will see below, the de- 
pendences in Figs. [6] and [7] obtained under the proposed 
theoretical approach, are qualitatively in good agreement 
with obtained experimental results. 



VI. EXPERIMENTAL PART: THERMAL 
EQUILIBRIUM FOR COUPLED ATOM-LIGHT 
STATES 



FIG. 7: (Color online) The reduced intensity of fluorescence I 
versus the square of the resonant Rabi frequency for 400-bar 
helium buffer gas at three different frequencies of the laser 
field. 



In this limit the maximal intensity weig ht is l( the 
(for S > 0) and AI approaches 



the 



r tanh (H \S\ /2k B T) « T. (42) 



On the contrary, for small Rabi frequencies (curves de- 
scribing nonequilibrium states in Figjf)]), from (1401) we 
have 



AI : 



S 2 



t&nh (h\S\/ 2k B T) <T, 



(43) 



which implies fulfillment of condition AI <C AI^ therm ^ 
[cf. ©]. 

In Fig J7] calculated dependencies of the total normal- 
ized intensity weight I versus the square of the resonant 
Rabi frequency (as a measure of the intensity) are pre- 
sented. Physically, they reflect the fact that the popula- 
tion of the atomic ground state \a) decreases, while that 
of the upper level (population abb) increases, with in- 



We now proceed to experimental work, in which we 
have investigated the rubidium atoms under an extraor- 
dinarily high buffer gas pressure driven by tunable laser 
radiation. The aim of our work presently is to reach ther- 
mal equilibrium of coupled atom-light states. Clearly, 
thermal equilibrium is a prerequisite for observation of 
a Bose-Einstein-like phase transition for polaritons. The 
large collisional broadening in our system interpolates be- 
tween the usual atomic gas phase and the solid or liquid 
phase condition. 



A. Experimental setup 

A scheme of our experimental setup used to investi- 
gate atom-light interaction at ultrahigh buffer gas pres- 
sures and intense laser radiation is shown in Fig 15] A 
stainless-steel high-pressure cell filled with atomic rubid- 
ium is used, which is connected by a valve to a buffer 
gas reservoir. Alternatively helium or argon buffer gas 
at pressures up to 230 bar is filled into the cell at room 
temperature. By heating the sealed cell up to 530 K, a 
buffer gas pressure of 500 bar for argon and a rubidium 
density of 10 16 cm -3 are reached. Tunable laser radia- 
tion derived from a titanium-sapphire laser tuned to the 
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FIG. 8: Experimental setup. 

rubidium D lines is focused on the cell to a waist size 
of 3/xm. The focus is placed directly behind the window 
to suppress propagation effects. To detect fluorescence 
light selectively from the focal region, where high laser 
intensities, up to 10 9 mW/cm 2 , for an optical power of 
300 mW are reached, both the incident beam and the 
outgoing fluorescence are spatially filtered with pinholes 
in a confocal geometry. 

B. Excited-state lifetime measurements 

In initial experiments, we have measured the lifetime 
of the rubidium 5P excited state in the presence of high 
buffer gas pressures. The lifetime of excited states is a 
major limitation of the coherence time in semiconduc- 
tor polariton BEC experiments, which is here typically 
close to 1 ps [1(. For the excited-state lifetime measure- 
ment in our rubidium high-pressure buffer gas system, 
we have chopped the exiting Ti:sapphire laser beam with 
an acousto-optical modulator (not shown in Fig(5]) and 
measured the fluorescence in the backward direction with 
a photomultiplier. Measurements of the 5P-5S fluores- 
cence decay time for various buffer gas pressures of ar- 
gon and helium are shown in FiglHl For low pressures 
the decay time is significantly longer than the natural 
5P lifetime of 27 ns [27j. We ascribe the slow decay of 
the fluorescence signal to energy pooling 33], where in 
excited- atom/excited-atom collisions two Rb atoms pool 
their internal energy to produce one ground-state atom 
and one in a higher excited state. A highly excited atom 
with a long lifetime can subsequently decay into the 5P 
state. Therefore energy pooling yields a reoccupation 
process of the 5P state that increases the observed SP- 
SS fluorescence decay time. For higher buffer gas pres- 
sures the observed decay time is reduced and tends to 
approach the expected 5P natural lifetime. We attribute 
this to a decline in the energy pooling process that in- 
volves highly excited rubidium states at high buffer gas 
pressures, which reduces the efficiency of the repopula- 
tion process at high pressure values. This interpretation 



FIG. 9: Observed fluorescence decay time of the 5P rubid- 
ium excited state versus buffer gas pressure. Data recorded 
for argon (squares) and helium (circles) gas together with ex- 
ponential decay curves fitted to the data for helium (dashed 
line) and for argon (solid line). 

of the measured decay times is supported by the observed 
strong blue fluorescence near 420 nm, which gives clear 
evidence of the presence of energy pooling, at low buffer 
gas pressures. Further, with increasing pressure the in- 
tensity of the blue fluorescence decreases considerably. 

To guide the eye, exponential decay curves have been 
fitted to the lifetime measurement data. The graphs ap- 
proach a decay time of 28,3 ns ± 1,9 ns for helium and 
29,6 ns ± 2,2 ns for argon at high pressure values. These 
values would coincide with the natural lifetime within 
our experimental accuracy. 

No evidence of a decline in the 5P lifetime, possibly 
to quenching, was observed within our present measure- 
ment accuracy up to the quoted pressure value of 200 bar, 
which was the maximum available buffer gas pressure in 
these data sets. In future, we plan to extend these life- 
time measurements toward higher buffer gas pressures. 

C. Measurements indicating thermal 
quasiequilibrium of atom-light states 

In subsequent measurements we have recorded spectra 
of the rubidium D lines at ultra-high buffer gas pres- 
sures. Typical fluorescence spectra recorded at 530 K 
and 400 bar helium buffer gas pressure for variable opti- 
cal powers are shown in FigfTOla). Since the fluorescence 
is measured in the experiment, the population of dressed 
states and population of excited atomic level abb [see (|39|) ] 
are of interest as well. 

At a moderate optical power (P = 25 mW) the usual 
pressure-broadened rubidium D\ and Di lines are visible. 
The resonant Rabi frequency f2 /27r is about 0.03 THz in 
this case and we deal with state ©, which is placed far 
from the region of the resonant atom-field interaction (see 
FigO) ■ For atom- field detuning <5/27r = — 11 THz and for 
collisional broadening 7^4,-/2^ ~ 3.6 THz the minimal 
thermalization time occurs for 500 bar argon buffer gas 
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FIG. 10: (a) Rubidium fluorescence signal and (b) results of 
extrapolation of measurement data to infinite laser intensity 
as a function of the incident frequency: both for helium buffer 
gas. Error bars were obtained from the extrapolation fits. 
Data were fitted with a Fermi-Dirac distribution, assuming 
thermal equilibrium of the two-level dressed-state system [cf. 
GDI- 



pressure and Ttherm ~ 37.4 ns in this case, which is of 
the order of spontaneous emission lifetime T spon t = 27 ns. 
Thus, we have not yet reached thermal equilibrium for 
such a state, represented by dashed (blue) curves in Figs. 
SandEl 

In particular, at the experimentally achieved maximal 
power P ~ 300 mW, the resonant Rabi frequency is 
fl /2n ~ 0.1 THz and the time of thermalization Ttherm 
approaches 10.8 ns (for 400-bar helium buffer gas pres- 
sure) and 3.37 ns (for 500-bar argon buffer gas pressure). 
The latter value is essentially shorter than the sponta- 
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FIG. 11: Fluorescence signal at 400-bar helium buffer gas as 
a function of input laser power for three different values of 
laser frequency. 



neous emission lifetime T spont . In this sense we can speak 
about achieving equilibrium for coupled atom-light states 
[state (2) in FigJ2] in the experiment (see FigJS]). 

Further diminishmcnt of the thermalization time is 
possible by increasing the incident optical power. No- 
tably, full thermalization of the coupled atom-light state 
is achieved for infinite laser power (cf. Fig|3J or full line 
saturation (see FigUTJ). 

For the purpose of comparison the expected fluores- 
cence signal at full saturation versus frequency is shown 
in Fig llOf b). It was derived at each frequency point by 
extrapolating the observed fluorescence to infinite laser 
intensity. Interestingly, this spectrum of the extrapolated 
fluorescence can be fitted well by a Fermi-Dirac function 
1(8) (solid line) [see (|H] . We interpret these measure- 
ments as evidence for approach of dressed-states ther- 
mal equilibrium at high drive laser powers (cf. FigjS]). 
Strictly speaking, only in this case does the Franck- 
Condon factor (overlapping of quasimolecular wave func- 
tions) not matter. 

Figure [TTJ shows the experimentally observed fluores- 
cence versus optical power for three different laser fre- 
quencies for the helium buffer gas data. Near a line center 
(crosses) the curve has already saturated at a relatively 
low power level to an intermediate fluorescence level. For 
significant detuning, the optical power at which satura- 
tion is achieved is higher, as visible for both blue (open 
circles) and red (filled dots) detunings. 

Interestingly, the saturation level of fluorescence is 
clearly different in both cases, with the blue (red) de- 
tuning leading to the largest (smallest) values. At 
full saturation, where spontaneous processes are negli- 
gible compared to stimulated processes, we expect the 
dressed-state populations to be in full thermal equilib- 
rium (cf. Figs. S] and |7||. 

During the course of the experiments, we found that 
the results of the spectral measurements were critically 
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dependent on the purity of the buffer gas used and pre- 
sumably, also, the decomposition of residual impurities. 
The measurements shown in FiglTUlwere carried out us- 
ing the buffer gas argon/helium 5.0 (supplier: Air Liq- 
uide). In contrast, we did not observe significant satu- 
ration of the fluorescence at high laser powers for buffer 
gas from a different supplier (argon/helium 5.0; supplier, 
Praxair). We attribute this to different residual impu- 
rity composition. In a pure inert buffer gas system this 
perturbation apparently does not occur at a significant 
rate, which reveals the remarkable elasticity of alkali gas 
atoms in collisions with inert gas atoms. To experimen- 
tally reach full thermal equilibrium also at very large 
laser detuning it is necessary to enlarge the coupling in 
the atom-light system. This can be achieved by use of a 
higher laser power or a cavity-based system. 



D. Toward thermalization of atom-light states with 
a continuum of modes 



• 200°C 




wavelength [nm] 

FIG. 12: Relative optical intensity measured in the forward 
direction after the cell versus frequency for different cell tem- 
peratures. The value of the incident laser wavelength (790 
nm) is indicated by the dashed verticalline. 



Up to now thermal equilibrium of a coupled atom- 
light system has been investigated for the case of a two- 
level dressed-state system. The observation of a Bosc- 
Einstcin-like phase transition from a thermal state to an 
ordered polariton state, however, requires a continuum of 
modes (see, e.g., [34]]-[36j]), as, for example, provided by 
the transverse modes of an optical cavity. Alternatively, 
a suitable optical waveguide could be used. To allow 
accumulation of polariton modes at the lowest available 
energy, it is clear that a process must exist that allows 
a continuous change of photon wavelengths toward lower 
energies, at least in the region around the lowest available 
photon energy of the optical resonator. 

In our rubidium high-pressure buffer gas system, we 
have tested for the presence of such frequency shifting 
processes. For this experiment, we spectrally analyzed 
the radiation detected at the cell output in the forward 
direction behind a second window of the cell (see Fig|S]). 
Note that in this forward-direction measurement, be- 
cause of the comparatively long optical path through the 
cell (5 cm), the confocal parameter used (~ 70 fim) is 
much shorter than the interaction length, so that the 
beam diameter increases considerably at the end of the 
cell and the optical intensity there is comparatively low. 
Due to the long interaction length, we also expect prop- 
agation effects to play a role. Typical spectra, mea- 
sured with a grating spectrometer, are shown in FigfTSl 
for various cell temperatures. The incident laser wave- 
length here was 790 nm, that is, between the rubidium 
D\ and the rubidium Di lines. At moderate tempera- 
tures (200°C) the observed fluorescence spectrum covers 
the spectral range from 750 nm to 950 nm, with clear 
dips near the rubidium D lines, which we attribute to 
the greater absorption there. At higher temperatures 
and therefore higher Rb densities, the cell becomes opti- 
cally thick over a large spectral range near the rubidium 
resonances, and in this wavelength range little transmit- 



ted optical power is observed. The spectra show a broad 
fluorescence peak with a maximum wavelength exceed- 
ing 900 nm. The fluorescence signal observed at the 
cell output has maximum spectral intensity on the red 
side of the rubidium spectrum; that is, Stokes scatter- 
ing is clearly stronger than anti-Stokes processes. We 
interpret the results of these measurements as evidence 
that processes exist which could support thermalization 
of photonic modes also in the case of a continuous mode 
spectrum, as required for experimental observation of a 
BEC-like phase transition of polaritons. Clearly, the ef- 
ficiency of these wavelength shifting processes, and also 
thermalization of polariton modes, needs to be studied 
in detail in future work. For the rubidium buffer gas 
system, a reasonable choice of wavelength for the lowest- 
cavity mode seems to be between 800 nm and 900 nm, 
where optical modes are populated by Stokes scattering. 
The presence of the phase transition would change the 
broad fluorescence signal in the forward direction to a 
sharp peak located at the frequency of the lowest mode. 

VII. CONCLUSIONS AND OUTLOOK 

In this paper have we studied the thermalization of 
coupled atom-light states under the influence of OCs. 
We found that thermal equilibrium is possible by con- 
trolling the resonant Rabi frequency and the atom-light 
detuning. We have shown that a nonvanishing macro- 
scopic polarization of the atomic medium occurs when 
the problem is described without the secular approxima- 
tion, which is typically used in the problem under dis- 
cussion. Experimentally we find evidence for a thermal- 
ization of the dressed atom-light states in an ultrahigh- 
pressure buffer gas environment. The observed intensity- 
dependent asymmetry of the spectra is interpreted as 
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partial thermal equilibrium of dressed atom-light states, 
when the driving-field detuning is chosen as \S\ ~ ksT/H. 
The thermalization process results in significant energy 
redistribution within the two sideband intensity compo- 
nents. We have observed Stokes scattering and charac- 
terized the lifetime of excitations in the presence of the 
buffer gas. 

In the future, it will be important to add a spatial 
confinement to allow dispersion of the coupled atom-light 
eigenstates with a low- frequency cutoff, suitable for BEC. 
This can be implemented using either a resonator or a 
waveguide structure. This also yields an enhancement of 
the field amplitude, which can allow full thermalization 
of coupled atom-light states in the high-pressure buffer 
gas system. When the rubidium density is increased, a 
strong coupling limit should be achievable in the buffer 
gas system, making polaritons relevant atom-light exci- 
tations. On the theoretical side it will be important to 
extend the treatment of coupled atom-light state (polari- 
ton) thermalization in the presence of OCs, as done in 



the present work. 

We conclude that the ultrahigh-pressure buffer gas ap- 
proach is a promising candidate for possible realization 
of a BEC-like phase transition of polaritons in an atomic 
physics system. 
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